Non adiabatic quantum search algorithms 
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We present two new continuous time quantum search algorithms similar to the adiabatic search 
algorithm, but now without an adiabatic evolution. We find that both algorithms work for a wide 
range of values of the parameters of the Hamiltonian, and one of them has, as an additional feature 
that, for values of time larger than a characteristic one, it will converge to a state which can be close 
to the searched state. 



I. INTRODUCTION 

Quantum computation has attracted the attention of 
researchers from several different areas This field 
of knowledge presents new scientific challenges to learn 
how to work with quantum properties to obtain more 
efficient algorithms. However, relatively few quantum al- 
gorithms were created; among them, Shor's and Grover's 
[2, ba] algorithms are the best known. Grover's search 
algorithm locates a marked item in an unsorted list of 
N elements in a number of steps proportional to y/N, 
instead of O(N) as in the classical case. It performs a 
unitary transformation of the initial quantum state so 
as to increase the likelihood that the marked state of 
interest will be measured at the output (amplification 
technique). It has been proved that there are neither 
quantum nor classical algorithms that can perform faster 
such an unstructured search This search algorithm 
has also a continuous time version []| that has been de- 
scribed as the analogue of the original Grover algorithm. 
From this continuous time version, and using the quan- 
tum adiabatic theorem^ adiabatic search algorithms have 
been developed 0, S, [a], that consist in guessing a time- 
dependent Hamiltonian whose dynamics evolves slowly 
enough so that it remains always near its instantaneous 
ground state. They solve the search problem in a time 
proportional to y/N /5, where S is a precision parameter 
that depends on the energy difference between the two 
lowest states. 

Another way to generate a continuous time quantum 
search algorithm [10j has been recently developed, that 
finds a discrete eigenstate of a given Hamiltonian Ho. 
This algorithm behaves like Grovers's, and explicitly 
shows that the search algorithm is essentially a resonance 
phenomenon between the initial and the searched states 

In this work we present two new continuous time 
search algorithms that are controlled by a time depen- 
dent Hamiltonian, similarly to the case of the quantum 
adiabatic search algorithm, but now the evolution is non 
adiabatic; then it is not necessary to impose slowness to 
the dynamics in order to preserve the system in the fun- 



damental state. These algorithms provide new insights 
to search algorithms, in particular: a connection between 
the resonant and the adiabatic search algorithms (for the 
first case), or the possibility to generate a new type of 
search algorithm in which one does not need to pick up 
a particular instant of time when the measure has to be 
performed (for the second case) , provided the parameters 
characterizing the Hamiltonian are conveniently chosen. 
In this second case, one reaches an asymptotic form for 
the searching probability, which is rather independent on 
the size of the database N, at the cost of increasing the 
energy resources in the Hamiltonian as N grows. 

The paper is organized as follows. In the next two 
sections we develop the two models of non adiabatic al- 
gorithms. In the last section we draw the conclusions of 
this work. 



II. NON ADIABATIC ALGORITHM I 

Consider N items in the database, each associ- 
ated with a vector in the complete orthonormal set 
{|n), n = 1,2, —,N} in a Hilbert space. Let us call \s) 
the unknown searched state that is associated with the 
marked item belonging to the previous group. We as- 
sume that the initial state is the symmetric normalized 
state 



|V>o 



i N 

^y>> 



and define the two Hamiltonians 



H = I-\i> ) W>o|, 



(1) 



(2) 



(3) 



where I is the identity matrix, their ground states being 
IV'o) and \s) respectively. The algorithm is built on the 
following time-dependent Hamiltonian 



H(t) = f(t)H + g(t)H s 



(4) 



2 



where f(t) and g(t) are time dependent functions that 
will be defined later. Notice that H s plays the equiva- 
lent role to 'marking' the searched state in Grover's algo- 
rithm. The goal of the search algorithm is to change 
IV'o) into | a) or some approximation there of, follow- 
ing the dynamics generated by the Schrodinger equa- 
tion. In this problem, we can restrict the analysis to 
the two-dimensional space spanned by \s) and \p) = 

N 

, }- \n). The wave function is then expressed 



as 



\^{t)) = a 8 {t)\s) + a p {t)\ P ) , 



(5) 



for some o s (t ) , a p (t) such that | a s (t ) | 2 + 1 a p (t ) | 2 = 1 with 
a s (0) = jijN and a p (0) = y/l -1/N. In the {\s}, \p)} 
basis we have the following matrix for the Hamiltonian 



H(t) = 



(N-l)f(t) 

-VF=T/(t) 



-VN=Tf(t) 

f(t) + Ng(t) 



The above matrix can be rewritten under the form 



H(t) = ~(f + g)I+~u,(t)n-a 



(6) 



(7) 



wherei^) = y/{f - gf + Afg/N, n = 
-^(-2s/N^lf/N,0,f - g - 2f/N) is an unitary 
vector and a stands for the Pauli matrices. Defining 
n = (sin 9, 0, cos 6), one can easily obtain the functions 
/ and g as a function of 9, with the following result: 



/(*) = - 

g(t) = 



N 



N 



Lu{t) sin 9 (8) 
w(t) cos(6> + 0) (9) 



with sin (3 = (2 - N)/N, cos/3 = 2y/N^l/N. 

The first term in Eq. ([7]) is proportional to the identity, 
and therefore amounts to a common (time dependent) 
phase that can be ignored if one only wants to evaluate 
probabilities. We will concentrate on the second term, 
which has eigenvalues E±(t) = ±|w(t) with correspond- 
ing time-dependent eigenvectors \E + ,t >= (cos |, sin |) 
and \E-,t >= (— sin |, cos |), with respect to the ba- 
sis {\s},p}}. In this form, it becomes evident that the 
evolution originated from the Hamiltonian amounts to a 
(time-dependent) rotation in the space spanned by the 
states {|s),p)} with the goal of maximizing the probabil- 
ity of the \s) state. 

The Schrodinger equation in this basis, in units such 
that H = l, becomes 



da s (t) 

dt 
da p (t) 

dt 



iH(t) 



a s (t) 
a p (t) 



(10) 



Now we take the following steps: first, we change to a 
new basis where this Hamiltonian is diagonal; second, we 



solve the Schrodinger equation in that basis, and third, 
we return to the original basis, where we are searching the 
state \s). The wave function, \ip(t)) can be expressed as 
a combination of the time-dependent eigenstates \E + ,t), 
\E-,t), 



- a+(t)exp(-i J E+{t)dt)\E+,t) 
o 

/, 

+a_(i)exp(-i J E-(t)dt)\E-,t). (11) 
o 

We have two expressions for the wave function, one in 
the {|s), \p)} basis, Eq. |(5]), and the other one in the 
{\E+,t), \E-,t)} basis, Eq. $U§. The relation between 
both basis can be expressed as a relation between its 
coefficients, that is 



a s (t) 
a p (t) 



U\t) 



a + (t) exp(-i J E + (t)dt) 
o 

t 

a-(t) exp(-i / (t)dt) 
o 



(12) 



where 



cos 77 — sin ■ 



(13) 



The Schrodinger equation in the new coordinates is 



da+ (t) 
di 

dt 



-M\t)U{t) d ^M{t) 



til 



a+(t) 
o-(*) 



(14) 



where M n (t) = exp(-i J E+{t)dt), M 22 (t) = 

o 

t 

exp(-i / E+(t)dt) and M 12 {t) = M 21 (t) = 0. From here, 
o 

one easily arrives to 

da+(t) 



dt 
da-(t) 
dt 



= -fl(t) a_(i), 

= n*(t) o+(t), 



where 



fl(t) ee -i^exp(z / u(t)dt), 



(15) 
(16) 

(17) 



with 6 = Alternatively, we can rewrite Eq. (fl7|l as 

m - ^( d -^nt)-g(t) d m 

Nuj 2 \ dt Jy ' yv 1 dt ) 
t 

(18) 



x exp(i / uj{t)dt); 
o 
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we shall take /(O) = 1 and g(0) = in this section, then 
the initial condition in the new coordinates are a+(0) = 
and <2-(0) = 1. Up to this point, our treatment of the 
problem is similar to that of the adiabatic algorithm. 
Now to proceed further we shall choose the function 8(t) 
(or, equivalently, the functions /(f) and g(t)) for the non 
adiabatic approach that has similarities with the resonant 
search algorithm. As seen in Eq. I|17p . if we choose 9(t) 
so as to cancel the time dependence of the modulus of fl, 
the system Eqs. lfl5lfT6|) will have an oscillatory solution 
between the amplitudes a+(t) and a_(i), with a period 
proportional to Vn, which can be identified with the 
Grover search time. This situation reminds us of the 
resonant search algorithm [13] , but now these amplitudes 
are not the amplitudes a s (t) and a p (t). 

In order to solve analytically the system Eqs. ([151 116] ) 
we shall impose the conditions: 



dg(t) 

dt 



/(*) - 9{t) 



dt 



e {at + jY 



uj(t) = \at + j\ . 



(19) 



(20) 



In these expressions £ is a coupling parameter between 
the states a + (t) and a_(t), a is the velocity parameter 
associated to the energy gap and 7 = 1 is dictated by our 
choice of /(0) and g(0). One can check that Eqs. (|19|20p 
are equivalent to imposing 



and 



T) = 1- 



-"0 



(26) 



The solutions of the Eqs. l(24|) are the parabolic cylinder 
function D ri (z) In our case the general solution is 



W(z) = AiD v (z)+AaD n (-z), 



(27) 



where the coefficients are determined by the initial con- 
ditions a_(0) = 1, and da ^°^ = 0. These coefficients 
are 



A x = 



Ao = 



D v -i(-z ) 



D v -i(zo)D v (—zq) + Dr,(zo)D v -i(—zo) ' 



D v -i(z ) 



D v -i(zo)D v (-z ) + D ri (z )D v ^ 1 {-z ) ' 



(28) 



(29) 



Finally, the amplitude a+(t) can be calculated using the 
above result for a_(f) and Eq. lfl6|) . 

Let us discuss in more detail the qualitative behavior of 
the results we have obtained so far. For large N and finite 
e in such a way that fio — -C 1 it follows from the 

Eqs. {MlEaEfilEaEHlESl) that \a-(t)\ ~ 1, \a+(t)\ ~ 0. 
On the other hand, for large N one can approximate 9 ~ 
2flot. Then using Eq. ([12]) . the following approximation 
for the probabilities of the searched and the orthogonal 
states are obtained 



P s (f)~ sin 2 (n t), 



(30) 



e = 2n (21) 

with Slo a constant, which we rewrite as flo = 
this way conditions (|19I20[) simply imply both a mixing 
angle 0(t) and a gap uj(t) which evolve linearly with time, 
and thus determine the time evolution of functions /(f) 
and g(t), see Eqs. (|8l9)l . 

We can decouple Eqs. fT5l fT6 l) to obtain a differential 
equations for a_(f) (and similarly for a+(f)) 



dt 2 



da- 
~df 



Qn 2 a_ 



(22) 



This equation can be solved in the same way as was 
done in [13| . The change of variable 



W — exp 



-(af 2 + 2 7 f) 



a_, 



leads to 



where 



z = \fa exp(— i7r/4)(f + j/a) 



(23) 



(24) 



(25) 



p p (t) ~ cos 2 (n t), 



(31) 



which are valid whenever f satisfies (at + 7) > 0. Note 
that the Eqs. ((30l l3Tj) are independent of the value of 
a if the previous conditions are verified. Then, if we let 
the system evolve during a time r = , and we mea- 
sure immediately after that, the probability to obtain the 
searched state is equal to one. In this case our method be- 
haves qualitatively like Grover's. The parameter e allows 
us for a faster search (relatively to the standard Grover's 
algorithms): one can even obtain a characteristic time 
t ~ 1. This speedup is allowed because the energy scale 
in the Hamiltonian, defined by the functions /(f) and 
g(t) is large enough (c.f. Eq. (22) in [12]), provided 
that Oo <§C 1. For concreteness, we will adopt the value 
c = 1. Additionally we have recovered our interpretation 
of the search algorithm as a quantum resonance between 
states ficl [TBI, fla. fl7l] : now the resonance is between the 
searched and the orthogonal states. 

The above result shows that the non adiabatic algo- 
rithm works correctly for a > (remember that 7 = 1). 
We have verified this situation for several values of N 
and a > using the exact equation in Fig. [T] The figure 
shows a periodic behavior with the Grover characteristic 
time and the correctness of the approximation in Eqs. 
([30] l3"T]l as N increases. 
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Figure 1: Probabilities as a function of time in units of the 
characteristics time r = n ^ . The the full line corresponds to 
the searched state P a , the dashed line to the orthogonal state 
Pp. In all cases a > 0, but the probabilities are independent of 
the value of a. The sizes of the searched set are: (a) N = 50, 
(b) N = 500, (c) N = 5000. 



When a = (see Eq. (|22|) ) a_ and a + can be easily 
obtained analytically 



o_(t) 



ex P(-*7^) 
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cos (/>(£) 
sin 0(i) 



a+(i) 



n exp(-i^t) 



sin 0(f) 



(32) 



(33) 



with (/>(t) = y (3Q + Oq£. From these expressions, and 

for large JV, f^o <C 1, |»+(£)l — 0) l a -(^)l — 1 and using 
the same arguments as before, it can be shown that the 
search algorithm also works in this case. For a < the 
behavior of the system is quite more complex. If, addi- 
tionally, | a | ~ (then both |zq| and \z\ go to infinity) 
using Eqs. (TJU E3 EH E2l S2l and the asymptotic prop- 
erty of the parabolic cylinder functions, it can be shown 
that P-(t) = |a_(t)| ~ 1 as before, then for large N and 



£7q < 1 Eqs. ([30] [31]) are again obtained, and the search 
algorithm continues to operate. For the case a < but 
finite we shall use another reasoning that could have also 
be used in the previous cases. Notice that the character- 
istic frequencies of the probability amplitudes (~ Qo) are, 

in general, very small compared with the time-dependent 

t 

characteristic frequency of O(t) i.e. \ J uj(t)dt, then the 

o 

stationary phase method can be used to integrate ap- 
proximately the differential equations Eqs. lfl5l fl6)) for 
a 7^ 0. We have used this method in the case a < 
and a finite, to obtain Eqs. f30l l3Tjl with the condition 
t < t c = — 2 The time t c is the 'close approach time', 
defined as the time when the derivative of the phase of 
Q(t) vanishes (see Eq. ((T7|l ) and at the same time the 

energy levels cross each other. Then the search algorithm 

t 

operates up to this time if f^o <C \ J ui(t)dt. Fig. [2] was 

o 

obtained using the exact results of this paper. It shows 
the probability P s for several values of a, and also that 
the approximation made in Eqs. l(30l [31]) remains valid 
for times t < t c . For times t > t c our previous argument 
cannot be applied, and the periodicity of the behavior 
is not clear because other frequencies are present. This 
figure establishes that, as \a\ is decreased, the close ap- 
proach time t c increases; in the limit, the algorithm works 
for all times. 

To close this section, let us remark that the previous 
results can also be obtained without the use of the mobile 
base in the following way: Let us substitute the expres- 
sions of f(t) and g{t) given by the Eqs. (|8|9|) in the 
equations Eq. (fTO]) . These equations, together with the 
normalization of the wave function and the conditions 
of maximization of the amplitude of the searched state, 
allow to find the time for which the probability of the 
search state is maximum. The maximization conditions 
= y |a s (*)| = 1 and, as a result, we obtain 



are 



da s {t) 



dt 



that they are satisfied for those values of t such that 



fir) 



N (err + 7) 



cos {2n T + (3) = 0. (34) 



From this equation it is deduced that, if ar + 7 7^ then 

nN , 



- le 



(35) 



with I = 1, 2, 3, .... As a result, we see that the algorithm 
works in an equivalent way to the Grover algorithm for all 
times, if a is positive, and only until the close approach 
time if a is negative. 



III. NON ADIABATIC ALGORITHM II 

In this section we introduce a new idea for the search- 
ing Hamiltonian which is based on a different choice of 
the functions fit) and g(t) and possesses the character- 
istic that one does not need to single out a given time in 
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Figure 2: Probability of the searched state P a as a function of 
time. N = 5000, in (a) a = -O.3I7/T, in (b) a = -O.IO7/T 
in (c) a = — 0.057/r. Note that the values t c change inversely 
with a. 



order to find the searched state with a high probability, 
provided the parameters of the Hamiltonian are chosen 
appropriately. 

Let us return to Eq. © and choose now f(t) = , 
then we have: 



With the above definitions, apart from a global phase 
which we will ignore, the Hamiltonian Eq. (|36| gives rise 
to the same evolution as the matrix 



H'(t) = 



-1 
-1 2b-2at 



(38) 



We will allow time to run from t = to arbitrarily large 
values (t — > 00). As we observe, the above Hamiltonian 
bears a close resemblance to the usual ones introduced 
in adiabatic quantum computations, in the sense that it 
has a time variation which is linear in time. However, in 
our case we will not start from the ground state of the 
Hamiltonian, and we will not intend either to force the 
system to be driven to its ground state for some finite 
time T by making use of the adiabatic theorem. 

The resulting evolution equations for a s (t) and a p (t) 
can be easily decoupled, leading to 



d a s 



„ „ s da s 
2i(b-at)—^ 
y ' dt 







(39) 



This equation has to be s upplemented by the initial 
conditions a s (0) = y/T/N , %fc{t = 0) = ia p (0) = 
iy/(N- Ty/N. With the substitution 



W = a s exp[-i(ia< 2 - bt)] 



(40) 



we arrive to the same equation as in Eq. (|24|) . but now 
2 = \j2aiit — b/a) and 77 = —jz. The solution to this 
equation can still be written in the form (|27|) . with co- 
efficients A\ and A 2 which have to be determined from 
the initial conditions. After some algebra, we arrive to 



Ai = 



1 



A,_i(-ab) + VNqoD^-za) 



N D n -i(z )D ri {~z a ) + D n {zQ)D n -\(-zo) 



(41) 



H(t) = y/N -II + 





-1 




(36) 



The first term in this expression is constant and pro- 
portional to the identity. One can, as done in the 
previous section, ignore it for the sake of solving the 
Schrodinger equation. Let us choose the (2,2) matrix 
element in Eq. ([361 so that it changes linearly with time. 
In this way, the resulting Hamiltonian (in the fixed ba- 
sis) mimics the evolution of the functions a+(t) and a_(t) 
obtained in Section 2. To be more precise, we choose 



N -2 



2(6 -at), 



(37) 



with a and b constants. Notice that f(t) and g(t) scale 
with N in the same way as in the previous section, there- 
fore the same discussion regarding the energy cost will 
apply. In this second model, the gap ener gy function 
takes the simple form co(t) — \J (at — b) 2 + 1. 



.4, 



D v -i(z ) 



Nq a D v (z ) 



N D, 7 _i(zo)A 7 (-.zo) + Di)(zo)D v -i(—zo) ' 



(42) 



exp(i37r/4) and 



20 



where qo = 

— byf^ exp(i7r/4). In order to give a result for the 

searched probability P s (t) — |a s (t)| 2 , we need to par- 
ticularize the values of TV, a and b. Using the asymptotic 
form for the parabolic cylinder functions [14]], one can 
obtain the following result for the limit lim^oo P s (t): 

p(a, b) = lim P s (t) = \A ie k ^ 4 + A 2 e-^ k ' i \ 2 1 (43) 
t — ^00 

with k = 

2 a 

Fig. [3] shows a contour plot of the limiting values of 
the probability p(a,b) with N = 100. As readily seen, it 
reveals a complicated pattern with bands of high proba- 
bility and low-probability valleys. These patterns depend 
quite weakly on TV and, in fact, it is possible to obtain 
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Figure 3: Contour plot for the limiting probability p(a, b) 
when N = 100. The horizontal axis corresponds to a, and 
the vertical axis to b. Brighter regions represent a higher 
probability, while dark regions indicate a probability close to 
zero. 
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Figure 4: Same as Fig. [3l for the asymptotic N — » oo limit. 



the limit N — > oo in Eqs. |42|) . Fig. [4] corresponds to this 
limit. As can be seen, the changes are moderate, showing 
that the asymptotic probability saturates for large values 
of N. 

In order to explicitly show the differences between our 
proposal and the adiabatic treatment, we have calculated 
the probability P s (t) for different parameters. The re- 
sults are shown in Fig. [H Given the structure seen in 
Fig. IU we choose a fixed b — 4.5, and plot the probability 
for three values of a: 1,5 and 20. The value of N we used 
is = 10 6 which, according to the previous discussion, is 
equivalent to taking N — > oo. 




3 6 9 12 15 

Time 



Figure 5: Evolution of the probability for the searched state 
P s as a function of time, for a fixed b = 4.5 and three values of 
the parameter a: 1 (solid line), 5 (thick dashes) and 20 (thin 
dashes) . 



In this figure, it is apparent that a transition occurs 
at the time t c = b/a, corresponding to the minimum dis- 
tance between the two eigenvalues of the Hamiltonian. 
The asymptotic behavior with time is clear from this fig- 
ure, although the final probability strongly depends on 
the particular choice of a and 6, as seen before. It is 
important to note that the corresponding Grover time 
would be of the order ViV ~ 10 3 , which represents a 
much larger time scale than the one showed in this fig- 
ure. 

Therefore, our proposal consists on implementing the 
Hamiltonian [6] but with the functions / and g as defined 
above. By an appropriate choice of the parameters a 
and b appearing in git), one can make the probability of 
the searched state to reach a value close to unity within a 
scale of time which is much shorter (for large N) than the 
corresponding Grover's time. At first sight, this appears 
to be in contradiction with the well-known result that 
Grover's algorithm is optimal for quantum searching 
However, it has been discussed that adiabatic search can 
be done at a shorter time (even with a time scale which is 
independent of N) at the expense of increasing the energy 
resources [l2l . [T3 |. This, in fact, seems to be the case 
within our proposal, if one remembers that, in order to 
obtain the Hamiltonian [36] using the resources defined in 
[6j the functions / and g will scale as V~N, thus increasing 
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the energy resources as N grows. 

IV. CONCLUSIONS 

We have developed two new continuous time quan- 
tum search algorithms using a time-dependent Hamil- 
tonians in a non adiabatic regime. Our approach dif- 
fers from to the usual (adiabatic) approach, when one 
starts from the initial ground state and tries to evolve 
slowly, making use of the adiabatic theorem to stay close 
to the instantaneous ground state. For the first case, 
the key of the algorithm is that the derivative of the 
amplitudes a_(t) and a+(t) have a fast time variation 
with a vanishing mean value over the characteristic time 

t i.e.i J — jfe — dt ~ 0, then starting from the ground 
o 

state, in the mobile basis {|-E+, t), \E- , i)}, the system 
remains near this ground state for all times. This algo- 
rithm behaves like the Grover algorithm for non negative 
values of the parameter a, independently of its particular 
value, for a large N and f2o 1- The optimal search 
time is proportional to yN, and the probability to find 
the searched state oscillates periodically. For a < the 
algorithm does not work for t > t c , with t c the close 
approach time. 



The second algorithm makes use of similar resources to 
build up a Hamiltonian that changes linearly with time. 
In our proposal, the initial and final states do not corre- 
spond to the ground states of the Hamiltonian, and the 
system is allowed to evolve up to arbitrarily large times, 
showing a convergence towards a final state after a fi- 
nite transition time. When the parameters are chosen 
appropriately, the asymptotic state can overlap with the 
searched state with high probability, and one does not 
need to pick up a special value of time to perform the 
measurement in order to obtain the desired result. 

Both algorithms can be used to perform a search within 
a time which can be made shorter than the standard 
Grover's time, at the expense of using also larger than 
standard energetic resources. 

These results open the possibility for the design of new 
quantum algorithms that perform a search on an unstruc- 
tured database (and possibly other algorithmic tasks) al- 
ternatively to the existing algorithms. 
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